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ABSTRACT 

Context. We present the results from the analysis of long Suzaku observations (149 ks and 122 ks for XIS and HXD, respectively) of 
the most X-ray luminous galaxy cluster, RX J1347.5-1 145, at z = 0.451. 

Aims. In order to understand the gas physics of a violent cluster merger, we study physical properties of the hot (~ 20 keV) gas clump 
in the south-east (SE) region discovered previously by the Sunyaev-Zel'dovich (SZ) effect observations. Using the hard X-ray data, a 
signature of non-thermal emission is also explored. 

Methods. We perform single as well as multi temperature fits to the Suzaku XIS spectra. Then the Suzaku XIS and HXD, and the 
Chandra ACIS-I data are fit jointly to examine the properties of the hot gas component in the SE region. Finally, we look for non- 
thermal emission in the Suzaku HXD data. 

Results. The single-temperature model fails to reproduce the 0.5-10 keV continuum emission and Fe-K lines measured by XIS 
simultaneously. The two-temperature model with a very hot component improves the fit, although the XIS data can only give a lower 
bound on the temperature of the hot component. We detect the hard X-ray emission in the Suzaku HXD data above the background 
in the 12-40 keV band at the 9cr level; however, the significance becomes marginal when the systematic error in the background 
estimation is included. With the joint analysis of the Suzaku and Chandra data, we determine the temperature of the hot gas in the 
SE region to be 25.3+^ (statistical; 90% CL) + _ 6 9 9 5 (systematic; 90% CL) keV, which is in an excellent agreement with the previous 
joint analysis of the SZ effect in radio and the Chandra X-ray data. This is the first time that the X-ray analysis alone gives a 
good measurement of the temperature of the hot component in the SE region, which is made possible by Suzaku's unprecedented 
sensitivity over the wide X-ray band. These results strongly indicate that RX J1347.5— 1 145 has undergone a recent, violent merger. 
The spectral analysis shows that the SE component is consistent with being thermal. We find the 3cr upper limit on the non-thermal 
flux, F < 8 x 10~ 12 ergs~'cnr 2 in the 12-60 keV band, which provides a limit on the inverse Compton scattering of relativistic 
electrons off the CMB photons. Combining this limit with a recent discovery of the radio mini halo in this cluster at 1.4 GHz, which 
measures the synchrotron radiation, we find a lower limit on the strength of the intracluster magnetic field B > 0.007 yuG. 

Key words, galaxies: clusters: individual: RX J1347.5-1 145 - galaxies: intergalactic medium - X-rays: galaxies: clusters - cosmol- 
ogy: observations 



1. Introduction 

RX J1347.5— 1 145 is one of the most extensively studied dis- 
tant clusters of galaxies, located at z = 0.451, which is also 
known as the most X-ray luminous cluster of galaxies on the 
sky. The bolome tric luminosity of RX J 1347. 5- 1145 is Lxboi ~ 
2 x 10 46 erg s" 1 (ISchindler et al.lll997l) . 

The previous multi-wavelength observations have shown 
that RX J 1347.5-1 145 has an unusually violent merger activity, 
which makes this cluster an ideal target for probing the intraclus- 
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ter gas physics and non-thermal phenomena associated with the 
cluster merger at high redshifts. 

The global temperature of the intracluster medium (ICM) 
of RX J1347.5-1145 is as high as 9-14 ke V, as indicated 
from observations with several X-ray satell i tes ( Schindler et alj 
IT9971 lEttoriet al.ll200lt lAllen et alj 12001 iGitti et all l2007bl) . 
This cluster has been classified as a "cooling-flow" cluster be- 
cause of its centrally peaked X-ray surface brightness profile as 
well as of its high mass accretion rate estimated at the center 
(ISchindler et al.ll 1997U Allen et aUl2002h . 

The most striking feature of RX J1347.5— 1 145 is the pres- 
ence of an extremely hot, ~ 20 keV, gas clump in the south- 
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east (SE) region. This component was discovered by radio 
observations of the Sunyaev-Zel'dovich (SZ) effect towards 
RX J1347.5-1 145 at 150 GHz as the prominent substructure 
of the SZ effect in the SE region, about 10" off of the center 
dKomatsu et al.l 120011) . iKitavama et all d2004l) have performed 
a detailed joint analysis of the SZ effect decrement data at 
150 GHz dKomatsuet all 120011) , increment data at 350 GHz 
dKomatsu et al.l Il999l) nind the Chandra ACIS-S3 X-ray data 
( Allen et al.ll2002lh and determined the temperature of the hot 
gas clump to be in excess of 20 keV, which is much higher than 
the average temperature of ambient gas. 

The coexistence of the cool and hot ICM in the central re- 
gion indicates a complex dynamical evolution of the system, 
such as a recent merger. In the optical band, the cluster center 
is dominated by two cD galaxies. The dark matter distribution 
has a lso been studied through the gravitational lensing effect 
(e.g. iFischer & Tvsonlll997l iSahu et al.lll998t IConen & Kneibl 
2002 :lAUen et al.ll2002H Bradac et al.ll2008t iMiranda et al.ll2 008; 
Halkola et alj 20081). An elonga ted, nearly bimodal distribution 



was found (BradaEetal]|2008), which also points to the clus- 
ter merger. Another line of evidence for the merger comes from 
a recent discovery of the radio mini halo at the cluster center 
dGitfietal] 12007a). which suggests the presence of relativistic 
particles powered by the merger as well as by th e cooling flow. 

The X-ray data from the ROSAT satellite dSchindler et alj 
1 1997b failed to identify the hot component in the SE region, as 
its X-ray sensitivity did not extend to such a high temperature. In 
fact, RX J1347.5— 1 145 had been considered as a relaxed, non- 
merging cluster until the SZ effect revealed the evidence for a 
violent merger. This example shows the importance of having 
sensitivity to high temperatures well in excess of 10 keV. 

Such a high temperature component seems common in merg- 
ing clusters. Another, perhaps most famous, example is the 
Bullet cluster (1E0 657-56), in which the h ot gas substructure 
has also been found ( Ma rkevitch et al.l2 002). The analysis of the 
temperature structure in merging clusters is a powerful tool for 
understanding the gas physics in extreme conditions. While the 
X-ray spectroscopy has been used widely for this sort of study, 
it is not very easy to precisely measure the temperature of very 
hot gas, as the majority of X-ray imaging spectroscopic obser- 
vations are limited to the energy band below 10 keV, and thus 
they cannot measure a spectral cut-off at kT characteristic of the 
thermal bremsstrahlung emission. 

To overcome this difficulty, Kitav ama et al.l (120041) used the 
fact that the SZ effect is sensitive to arbitrary high temperature 
gas, and performed a joint analysis of the SZ effect and X-ray 
images. They have determined the temperature of the SE excess 
component in RX J1347.5-1145 to be 28 + 7 keV, while the 
Chandra X-ray spectroscopy, whose sensitivity degrades signif- 
icantly beyond 7 keV, gave the lower limit of 22 keV. This is one 
of the hot test gas known t o date. Based on the nu merical simula- 
tions by Takizawa ( 1999), Kitava ma et al.l d2004) suggested that 
RX J 1347.5-1 145 has undergone a recent merger at the collision 
speed of 4500 kms _I , which is very simil ar to what has been 
found more recently from the B u llet cluster dMilosavlievic et al.l 
20071: IS^ringel & Farrarl 120071: iMastropietro & Burkertl 120071: 
Nusserll2008l) This result is indicative of RX J1347.5-1 145 be- 
ing essentially the Bull et cluster, excep t for th e viewing angle. 

The conclusion of Kitavama et al.l d2004l) relies on the SZ 
data from high-resolution mapping observation s of the SZ ef- 
fect ( t he hi ghest resolution ever achieved to date; Komats u et al.l 
[T999ll200lT) . However, precision of such SZ observations is still 
limited due to technological reasons; thus, high-precision X-ray 
observations whose sensitivity goes well beyond 10 keV are re- 
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Fig.l. Suzaku XIS-0 image of RX J1347.5-1145 in the 0.5- 
10 keV band. The spectral integration regions for the cluster and 
the background are indicated with the solid and dashed circles, 
respectively. The two corners of the CCD chip illuminated by 
55 Fe calibration sources are excluded from the image. 



quired to obtain reliable estimates of physical properties such as 
the temperature and merging velocity. 

The unpreceden ted sensitivity of the Suzaku satellite 
(Mits uda et al.ll2QQ~7h to hard X-ray emission well over 10 keV 
offers a great opportunity for studying violent merger events. 
In this paper we present an analysis of the temperature struc- 
ture and non-thermal, high-energy component in the distant clus- 
ter RX J1347.5— 1 145 with th e on-board instrum ents, the X-ray 



Imaging Spectrometer (XIS; Kovama et al. 2007 ), and the Hard 
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X-ray Detector (HXD; iTakahashi et al.l 120071; 
2007). These instruments provide high-sensitivity spectroscopic 
observations in a wide energy band up to several tens of keV. 
Combined with the Chandra data, we determine the temperature 
of the SE clump only from the X-ray spectroscopy, and discuss 
its physical properties. We also constrain, for the first time, the 
non-thermal hard X-ray emission and estimate the magnetic field 
strength in RX J1347.5-1 145. 

Throughout this paper we adopt a cosmological model with 
the matter density Q.m - 0.27, the cosmological constant £2a = 
0.73, and the Hubble constant Hq = 72 km s " 1 Mpc' 1 , which are 
consistent with the WMAP 5-year results dDunklev et al.| [2008: 
IKomatsu et~aT]|2008l) . At the cluster redshift (z = 0.451), 1" cor- 
responds to 5.74 kpc. Unless otherwise specified, quoted errors 
indicate the 90% confidence range. 



2. Observation and data reduction 

Long observations of RX J 1347. 5-1 145 were carried out twice 
in the Suzaku AO-1 period. During the observations, XIS and 
HXD were operated in their normal modes. To prioritize the 
HXD effective area, the target was observed at the HXD nom- 
inal aim-point (i.e., 3'. 5 off-axis relative to the XIS) . The total 
exposure time after data filtering is 149 ks/122 ks for XIS/HXD. 
The summary of the observations is given in Table Q] 

The XIS consists of four X-ray CCD cameras: three front- 
illuminated CCDs (XIS-0, -2, -3: 0.4-12 keV) and one back- 
illuminated CCD (XIS-1: 0.2-12 keV), and covers a field of 
view of 18' x 18'. The "Spaced Charge Injection" option was 
not applied. The energy resolution was ~ 160 eV at 6 keV 
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Target 




Sequence No. 


Date 


Coordinates" 


Exposure [s] b 










RA [deg] Dec [deg] 


XIS PIN 


RXJ1347.5- 


1145 


801013010 


2006-06-30 


206.8560 -11.8093 


69661 56698 


RXI1347.5- 


1145 


801013020 


2006-07-15 


206.8558 -11.8095 


79126 64922 



a Pointing coordinates in J2000. 

b Net exposure time after data filtering. 



nificantly affect our spectral analysis: if we compare the above 
arf files with those produced assuming a point-like source or an 
exposure-corrected Chandra image as an input source distribu- 
tion, these three kinds of arfs agree within 3% and give con- 
sistent spectral parameters within the present calibration uncer- 
tainty of the Suzaku telescope (see also AppendixlAl. 

The hard X-ray spectra were accumulated from all units of 
the PIN diodes. The PIN detector background was subtracted 
by using the Non-X-ray Background (NXB) file£] provided by 
the HXD instrument team. Since precise modeling of the back- 
ground is crucial for measuring the hard X-ray flux, we describe 
the method of background subtraction in $4] In the spectral fit- 
ting, we use the PIN response function^ that is appropriate for 
the epoch of observations. 



3. XIS analysis: 0.5-10 keV 

In order to examine the global temperature structure of 
RX J 1347.5-1 145, we fit the APEC thermal emission model to 
the observed XIS spectra in the 0.5-10 keV band. 

The angular size of the X-ray emission is fairly compact 
and nearly 90% of the total emission comes from the central 
r < 2' region, which is comparable to the spati al resolution of 
the S uzaku XRT (Half Power Diameter ~ 2'; Serlemitsos et al.l 
120071) . We use the extraction radius of 5'; thus, more than 97% 
of the photons from RX J1347.5-1 145 are collected. 

3.1. Fitting with the single-temperature model 

The XIS spectra with high photon statistics allow us to mea- 
sure the global temperature in two different ways: (1) the APEC 
model fit to the 0.5-10 keV spectra, and (2) the use of the inten- 
sity ratio of two Fe emission lines, namely (He-like Fe Ka)/(H- 
like Fe Ka), as an indicator of the "ionization" temperature. By 
comparing these two, we examine whether the emission can be 
modeled by the single-temperature gas, or the data require more 
complex temperature structure. 

For the analysis (1) (APEC model), we fit the APEC model 
to the observed 0.5-10 keV XIS spectra, where the redshift and 
the Galactic hydrogen colu mn density were fixed at z = 0.45 1 
and 7V H = 4.85 x 10 20 cm" 2 dDickev & Lockmanll 19901). respec- 
tively. The metal abundances table in lAnders & Grevessd (fl989) 
was assumed. The spectral bins in the 1.7-1.9 keV range were 
excluded due to the large calibration errors in the Si-edge struc- 
ture. In order to take into account the uncertainty in the thickness 
of the contaminating material, the absorption column densities 
of oxygen and carbon (Nc and No) were included in the spectral 
model with the ratio fixed at Nq/No = 6. In addition, there is 
a potential uncertainty in the absolute energy scale on the order 



(FWHM) at the time of the observations. As shown in Fig. [TJ 
the X-ray emission from the cluster is clearly detected. The X- 
ray peak position (13:47:30.3, -11:45:00.9) is consistent with 
that obtained from the Chandra (13:47:30.6,-1 1:45:09.3) within 
the acc uracy of the attitude d etermination of Suzaku (typically 
< 20"; lUchivama et aT]l2008l) . 

The non-imaging HXD covers a wide bandpass of 10- 
600 keV with PIN diodes (10-60 keV) and GSO scintillator (40- 
600 keV). We shall use the PIN data in this paper. The narrow 
field of view of ~ 30' x 30' (FWHM) and the low background 
level of the PIN spectral data enable us to study the hard X- 
ray emission from the cluster up to several tens of keV, ideal 
for studying a violent merger event. The bias voltage from one 
of four high-voltage units was reduced from the nominal value 
of 500 V to 400 V. As a result the total effective area of the 
HXD/PIN decreased by about a few % above 20 keV. This effect 
is incorporated in the detector response functions. 

We use the cleaned event files created by the pipeline 
processing version 2.0, and perform the data analysis using 
HEASOFT version 6.3.1. The XIS data was filtered according to 
the following criteria: the Earth elevation angle > 10°, the day- 
Earth elevation angle > 20°, and the satellite outside the South 
Atlantic Anomaly (SAA). The HXD data was filtered according 
to the following criteria: the Earth elevation angle > 5°, the geo- 
magnetic cut-off rigidity (COR) > 6 GV, and the satellite outside 
the SAA. 

The XIS spectra were extracted from a circular region within 
a radius of 5' that is centered at the X-ray peak position (see 
Fig- ID- The background spectra were accumulated from the cir- 
cular region within a 3' radius which is almost free from the 
cluster emission. Since the background contributes 10-15% to 
the total spectra, the systematic error due to the background un- 
certainty is small. To confirm this, we have checked the position 
dependence of the background spectra by analyzing the blank- 
sky data obtained during the ~ 90 ks North Ecliptic Pole obser- 
vation on 2006-02-10, and found that it is smaller than 7% for 
all of the four sensors of XIS. Therefore, the background does 
not affect our spectral analysis within the statistical errors. 

The energy response files were generated by using the 
FT00LS task, xisrmfgen. In order to take in to account the vi- 
gnetti ng effect of the X-ray telescopes (XRT; ISerlemitsos et alj 
2007) and a decrease of the low-energy efficiency due to the 
contaminating material on the optical blocking filter of the 
XIS, the auxiliary r esponse files (arfs) were calculated by us- 
ing xissimarfgen (Ish isaki et alj|2007l) . Here, we assume that 
the X-ray surface brightness distribution of the cluster is rep- 
resented by the double yS-model: our reanalysis of the Chandra 
ACIS-I image yields the two core radii given by r c \ = 24 kpc and 
r C 2 = 94 kpc, the same slope parameter, /3 = 0.63, for both com- 
ponents, and a surface brightness ratio of S i /S 2 = 8.1 at the cen- 
ter (see Fig. [9}}). Since the spatial distribution of RX J1347.5- 
1145 is compact compared to the typical size of the Suzaku' s 
Point Spread Function (PSF), its detailed feature does not sig- 



1 http://www.astro. isas.jaxa.jp/suzaku/analysis/hxd/pinnxb/pinnxb_ver2. 0/ 

2 ae_hxd_pinhxnome2_20070914.rsp 
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Table 2. Single-temperature APEC model parameters fit to the XIS spectra taken by four sensors (XIS-0, XIS-1, XIS-2, XIS-3), 
and the result from the simultaneous fit to all of the sensors (XIS-0, 1, 2, 3). 



Sensor 


kT [keV] 




Z [solar] 




K* 




x z /dof 


XIS-0 


13.31(12.76- 


13.86) 


0.35(0.28 


- 0.40) 


1.37(1.37- 


-1.39)xl0-- 


352.1/298 


XIS-1 


12.23(11.54- 


12.83) 


0.36(0.30 


-0.41) 


1.41(1.39- 


-1.44)xl0~ 2 


345.0/298 


XIS-2 


12.87(12.40- 


13.32) 


0.34(0.29 


-0.38) 


1.40(1.38- 


-1.42)xl0- 2 


316.5/298 


XIS-3 


12.76(12.28 - 


13.24) 


0.30(0.25 


-0.35) 


1.45(1.43- 


-1.47)xl0- 2 


300.2/298 


XIS-0, 1,2,3 


12.86(12.61 - 


12.94) 


0.33(0.31 


-0.36) 


1.37(1.36- 


-1.39)xl0- 2 


1320.4/1198 



a The APEC normalization, K = / n e n H dV/(47TD 2 A (l + zf) [10 14 errT 5 ]. 
Table 3. The instrumental calibration parameters from the APEC model fit simultaneously to the XIS-0, 1, 2, and 3. 



Sensor 


Relative normalization 


N, [10 18 cm- 2 ] 


Gain offset [eV] 


XIS-0 


l(Fix) 




0.012(< 0.116) 


-12.0 (-20.0—6.0) 


XIS-1 


1.035(1.023- 


1.047) 


0.093(0.014-0.184) 


-16.0 (-22.5—10.2) 


XIS-2 


1.019(1.007- 


1.030) 


0.01 1(< 0.116) 


-9.8 (-16.5— 3.1) 


XIS-3 


1.049(1.038- 


1.060) 


0.000(< 0.100) 


-20.1 (-26.4 — 14.2) 



(a) (b) 

E ' ' ' ' ' ' \ ° f ' ' ' ' ' He-like Fe Ka' 




Energy [keV] Energy [keV] 

Fig. 2. (a) Observed XIS spectra (r < 5')- The spectra taken by four sensors, XIS-0 (black), 1 (red), 2 (green), and 3 (blue), are shown 
separately. The solid lines in the upper panel show the best-fitting single-temperature APEC thermal plasma model simultaneously 
fit to all of the sensors, convolved with the telescope and the detector response functions. In the bottom panel the residuals of 
the fit in units of the number of standard deviations are shown, (b) Blow-up of the XIS-0 spectrum in the 4-5.5 keV band. The 
single-temperature model fails to fit the He-like Fe Ka line. 



„ W (b) 

o r ' ' ' ' i ' ' ' ' i ' ' ' ' : a r— 




Energy [keV] 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

°5 10 15 20 

kT [keV] 



Fig. 3. (a) Observed XIS spectra in the 4-5.5 keV band. The spectra taken by four sensors, XIS-0 (black), 1 (red), 2 (green), and 
3 (blue), are shown separately. The solid lines in the upper panels are the best-fitting continuum plus two Gaussian models. In the 
bottom panel the residuals of the fit in units of the number of standard deviations are shown, (b) Line ratio (He-like Fe Ko-)/(H-like 
Fe Ka) as a function of temperature, calculated from the APEC model. The horizontal dashed line and the dotted lines show the 
best-fitting line ratio and the 90% confidence interval, respectively. 
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of 10-20 eV, as indicated from the analysis of the 55 Fe calibra- 
tion source spectra; thus, the gain offset for each XIS sensor is 
included in the fitting model. 

Since all of the four sensors yield consistent APEC param- 
eters (see Table |2|, we shall show the results from the simulta- 
neous fit to all of the sensors (XIS-0, XIS-1, XIS-2, and XIS-3) 
hereafter. When performing the simultaneous fit, the normaliza- 
tions of XIS-1, XIS-2, and XIS-3 relative to that of XIS-0 were 
also treated as free parameters. We find the temperature and the 
metal abundance of kT = 12.86^ keV and Z = 0.33+° ^ so- 
lar, respectively. The unabsorbed flux and luminosity in the 0.5- 
10 keV band are ^x,o.5-io keV = 1-3 x 10~ n ergs -1 cm -2 , and 
£x,o.5-iokeV = 8.7 x 10 45 ergs" 1 , respectively. The estimated 
bolometric X-ray luminosity is Lx.boi = 1-37 x 10 46 ergs -1 . In 
Table [3] we list the the instrumental calibration parameters. In 
the following analysis we shall use the best-fitting values for the 
relative normalization factors, the carbon column densities, and 
the gain offsets. 

Fig. [2] shows the observed XIS spectra. Fig. |2j3 is a blow- 
up of the region around the redshifted He-like Fe Ka and H- 
like Fe Ka lines at about 4.6 keV and 4.8 keV, respectively. We 
find clear detections of both lines. Also shown is the best-fitting 
single-temperature APEC thermal plasma model. The single- 
temperature model fails to fit the He-like line. The reduced^ 2 
from the single-temperature fit is x 2 = 1320.4/1198; thus, the 
single-temperature model is rejected at the 99.3% C.L. 

We then attempt to fit the XIS spectra with the NEI (Non- 
equilibrium ionization collisional plasma) model. We found that 
the NEI model does not improve the fit to the He-like Fe Ka line 
and the quality of the fit in the whole energy range is similar to 
the 12.9 keV APEC model: the resulting^ 2 from the NEI model 
fit is 1340 for 1208 degrees of freedom. Moreover, the best-fit 
ionization timescale of > 10 12 s cm 3 is a ctually long e nough for 
the system to reach the equilibrium state (Masai 1994), and thus 
the assumption that the system is in a non-equilibrium state can- 
not be justified. From these results we conclude that the single- 
temperature NEI model fails to describe the observed XIS spec- 
tra. 

For the analysis (2) (line ratio), we model the observed XIS 
spectra in the 4-5.5 keV band by the sum of the following com- 
ponents: (i) the APEC model with the metal abundance reset to 
for the continuum emission, (ii) two Gaussian functions for the 
two major Fe lines at 4.6 and 4.8 keV, and (iii) another Gaussian 
line at 5.4 keV for the blend of the He-like Ni-Ka and He-like 
Fe-K/3 lines. In Fig. [3^ we show the results of the fit. This model 
gives an acceptable fit to the data with^ 2 /dof = 185/163, or the 
probability of finding higher^- 2 values than observed is 1 1.4%. 

This model provides the line ratio, (He-like Fe Ka)/(H-like 
Fe Ka), of 1.05(0.92 - 1.34). Fig. [3p shows the temperature 
dependence of the line ratio predicted by the APEC model, 
folded with the XRT+XIS response functions. From this we find 
kT = 10.4(9.1 - 11.4) keV, which is significantly different from 
that found in the analysis (1), kT = 12.86(12.61 - 12.94) keV, 
suggesting that the ICM cannot be explained by the single- 
temperature model. Therefore, the ICM is more likely to be in 
the form of multi-temperature plasma, which we shall explore in 
§E2] 

Since the iron lines are emitted mostly from low-temperature 
gas (see j j3.2| >, the line energies and widths are useful probes of 
(i) the gas bulk motion and (ii) the turbulence in the cluster core. 
For (i), the observed centroid energies of the iron lines are con- 
sistent with those expected from the 10 keV APEC model (6.69 
and 6.97 keV in the rest frame) and the cluster redshift, to within 
the 90% statistical error of 10 eV. As shown in Table [3] there 



is the gain offset of 15 eV when averaged over the four XIS 
sensors, we simply assign the 90% systematic error of 15 eV 
to the XIS energy scale. Then by adding the statistical and sys- 
tematic errors in quadrature, the 90% upper limit on the line- 
of-sight bulk velocity relative to z — 0.451 is estimated as 
AV < 1200 kms -1 . (ii) We find no significant line broadening: 
the 90% upper limit on the Gaussian width of the He-like Ka 
line is cr measuret j < 35 eV. Since the measured Gaussian width is 
the sum of the thermal Doppler broadening, <x (<x = 4 eV for 
kT = 10 keV), and the broadening due to a turbulent velocity, 
o- turb , i.e., cr 2 ^^ = o- 2 + cr 2 mb , an upper limit on the turbulent 
velocity dispersion is obtained as cr tllrb < 1500 km s -1 . Therefore 
we did not find from the analysis of the iron lines any significant 
bulk velocity and the turbulent velocity, which is consistent with 
the view that the low-temperature gas in the cluster core is rel- 
atively relaxed. The derived upper limits are also com parable to 
those obtained for the core of the Centaurus cluster (lOta et all 
2007). If the velocity structure on a small (< 2') spatial scale ex- 
ists, the signals may be diluted due to the wide PSF of Suzaku. 
However, it is not easy to study the detailed spatial structure in 
RX J 1347.5-1 145 because the angular size is small. 

3.2. Fitting with the two-temperature model 

Given that the single-temperature model fails to describe the ob- 
served XIS spectra, we explore the two-temperature model con- 
sisting of two APEC models with different temperatures, as the 
simplest case of the multi-temperature plasma. The results of the 
fit to the two-temperature model are shown in Fig.|4]and Table|4] 
where the metal abundances of the two components are assumed 
to be equal. 

From the fit to the 0.5-10 keV spectra, we find that the 
low-temperature component has the temperature of kT\ 0W = 
9.V°{] ke\fl wher eas we only find a lower bound on the tem- 
perature of the high-temperature component as £7hi g h > 33 keV. 
This two-temperature model yields an acceptable goodness-of- 
fit: the reduced X 1 is x 1 = 1295.0/1207, or the probability of 
finding higher^- 2 values than observed is 3.9%. In comparison 
to the single-temperature APEC model, the fit has improved at 
more than the 99.99% confidence level according to the F-test 
(Ax 2 = 25 for two additional parameters). Restricting the en- 
ergy range to 4.0-5.5 keV does not change the parameters sig- 
nificantly (see the lower row of Table |4|. The reduced x 2 1S 
X 2 = 185.2/159, or the probability of finding higher^ 2 values 
than observed is 7.6%. 

Since the observed Fe Ka lines should be totally dominated 
by the low-temperature gas (see Fig. |4j)), one may compare 
£7k>w = 9.7^j ] keV and the temperature derived from the line ra- 
tio, kT = 1 0.4(9 . 1 - 1 1 .4) ke V (see § |3Tj, directly. These two es- 
timates are in an excellent agreement, which gives a support for 
the two-temperature interpretation of the ICM of RX J 1347.5- 
1145. 

Although the XIS spectra give the lower bound on the tem- 
perature of the hotter component (see Table 0), they fail to give 
the upper bound. We find that this is due to the two-temperature 
model being too simple to be realistic. More detailed model- 
ing on the hotter component, including a joint analysis with the 
Suzaku HXD data and the spatially-resolved Chandra data, will 
be given in § [5] 



3 The temperature of the low-temperature component is similar to the 
temperature determined from the spectra taken by the GIS detectors on 
board ASCA with the single-temperature model, kTAscA = 9.3*Ji keV 
<SchindleretaIlll997h . 
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Table 4. Two-temperature APEC model parameters fit to the XIS spectra in the 0.5-10 keV and 4.0-5.0 keV bands. 



Energy band 


;' 


kT t [keV] 


Z, [solar] 






V/dof 


0.5-10 keV 


1 


9.7(8.6-10.4) 


0.33(0.31-0.36) 


1.07(0.88 


- 1.13) x 10~ 2 


1295.0/1207 




2 


> 32.5 


= Zi 


3.73(3.01 


- 5.39) x 10~ 3 




4.0-5.5 keV 


1 


9.8(3.6-11.6) 


0.35(0.28-0.52) 


9.2(2.9 - 


14.1) x 10~ J 


185.2/159 




2 


34.4(> 22.2) 


= z, 


4.67(0.26 


- 1.07) x 10~ 3 






- (b) 





Energy [keV] 



4.5 5 
Energy (keV) 



5.5 



Fig. 4. (a) Two-temperature APEC model fit to the observed XIS spectra. The spectra taken by four sensors, XIS-0 (black), 1 (red), 
2 (green), and 3 (blue), are shown separately. The solid lines in the upper panel show the sum of two APEC models with different 
temperatures, while the dashed lines show the low-temperature (upper curves) and high-temperature (lower curves) components 
separately, convolved with the telescope and the detector response functions. In the bottom panel the residuals of the fit in units of 
the number of standard deviations are shown, (b) Blow-up of the XIS spectra in the 4-5.5 keV band. The two-temperature model 
fits both the He-like and H-like Fe Ka lines adequately. 



4. HXD/PIN analysis: 12-60 keV 

For the PIN diodes of HXD working in the energy band of 12- 
60 keV, the background is composed of the following three com- 
ponents: (i) the non-X-ray background (NXB), (ii) Cosmic X- 
ray Background (CXB), and (iii) bright point-like sources within 
the detector field of views. Of these, the NXB dominates the 
observed flux particularly for faint hard X-ray sources such as 
clusters of galaxies. 

Since the NXB has a significant time-variability, the repro- 
ducibility of the NXB is the most important factor in the mea- 
surements of the hard X-ray flux from clusters of galaxies. In 
§ 14.11 the systematic error due to NXB is carefully examined. 
In § I4.2l the CXB model is calculated. In §[43] the point source 
fluxes are estimated from the XMM-Newton observations of the 
same field. 

4. 1 . Non-X-ray background 

The latest report from the instrument tearrQ shows that the re- 
producibility of the PIN NXB model for the data that have been 
processed with the version 2 pipeline was typically 3% (lcr). 
This estimate is based upon systematic comparisons between the 
NXB model and the PIN data during the Earth occultation in a 
trend archive. We check below the validity of assigning the 3% 
systematic error to our PIN analysis. 

In Fig. [5] we show comparisons between the observed PIN 
data (hereafter "Data") and the NXB model in two observing pe- 
riods, (a) Seq: 801013010, and (b) Seq: 801013020. In Table 
we summarize the count rates in the 12^-0 keV and 40-60 keV 
bands. The excess signal has been detected in Data-NXB at en- 
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ergies below 30 keV during both observing periods. The dif- 
ference in Data-NXB between two observing periods is well 
below the NXB itself: (Data - NXB) 2 - (Data - NXB)i ~ 
0.02(NXBi + NXB 2 )/2. However, the difference is marginally 
(~ 2<t) inconsistent with zero compared to its statistical error: 
(Data - NXB) 2 - (Data - NXB)i = 0.0088 + 0.0043, which 
could be indicative of some residual time variability. 

In general, the NXB intensity changes depending pri- 
marily on the satelli te's passage of the SAA and the COR 
dKokubun et al.ll2007l) . To study these effects, we bin the PIN 
spectra under the following different conditions: elapsed time 
after the passage of the SAA, TSAA% 6000 s, and three dif- 
ferent COR [GV] ranges, 6 < COR < 10, 10 < COR < 12, 
12 < COR < 14. This analysis shows that there is a trend 
to obtain a higher count rate for TSAA< 6000 s compared to 
TSAA> 6000 s, and a higher count rate for a smaller COR. 
However, the differences in the count rates are below +3% of 
the NXB intensity in each case. 

In order to study the reproducibility of the NXB further, a 
comparison was made for Data and NXB during the periods of 
the Earth occultation. The Earth is known to be dark in hard 
X-rays and the Earth occultation data can be regarded as NXB 
for the PIN observations. As a result, Data and NXB during the 
Earth occultation show a good agreement, and Data-NXB is 
consistent with zero. From these studies we conclude that the 
systematic error of the NXB model is fairly small (< 1 %). 

In summary, the accuracy of the PIN NXB model is consis- 
tent with that found in the latest report from the instrument team. 
Therefore, we shall adopt the instrument team's estimate of the 
systematic error, 3%, and propagate it through our spectral anal- 
ysis of the HXD/PIN data. 
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Fig. 5. Observed PIN spectra during two observing periods of RX J1347.5-1 145. (a) Seq: 801013010, and (b) Seq: 801013020. In 
each panel Data, NXB, and (Data-NXB) are shown with the black crosses, red crosses, and green circles, respectively. 
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Fig. 6. Observed light curves of the PIN 12^-0 keV count rates during two observing periods, (a) Seq: 801013010, and (b) Seq: 
80101 3020. In the upper panels the Data and NXB rates are shown with the black and red crosses, respectively. In the bottom panels 
(Data-NXB) is shown with the green circles. The error bars show the lcr statistical errors. 



Table 5. PIN count rates measured in the 12^-0 keV and 40-60 keV bands. The quoted errors are lcr statistical errors. 



Seq 


Energy range 


Exposure 


Data 




NXB 




(Data-NXB) 


(Data-NXB )/NXB 






[s] 


[cs- 1 ] 




[cs- 1 ] 




[cs- 1 ] 




% 




801013010 


12^10 keV 


56698 


0.5180 


± 0.0030 


0.4754 ± 


0.0009 


0.0426 ± 


0.0031 


8.96 ± 


0.67 


801013020 


12-40 keV 


64922 


0.5223 


± 0.0028 


0.4709 ± 


0.0008 


0.0514 ± 


0.0030 


10.92 


+ 0.63 


801013010+801013020 


12^10 keV 


121620 


0.5203 


±0.0021 


0.4730 ± 


0.0006 


0.0473 ± 


0.0022 


10.00 


±0.46 


801013010 


40-60 keV 


56698 


0.0497 


± 0.0009 


0.0499 ± 


0.0003 


-0.0002 


± 0.0010 


-0.44 


± 1.96 


801013020 


40-60 keV 


64922 


0.0497 


± 0.0009 


0.0495 ± 


0.0003 


0.0002 ± 


0.0009 


0.43 + 


1.85 


801013010+801013020 


40-60 keV 


121620 


0.0497 


± 0.0006 


0.0497 ± 


0.0002 


0.0000 ± 


0.0007 


0.03 ± 


1.35 



4.2. Cosmic X-ray Background 

We calculate the CXB spectrum in the PIN band assuming a 
power-law model with an exponential c ut-off at 40 k eV as pre- 
viously reported from the HEAO-1 A2 dBoldtll 1987h . Since the 
PIN detector response assumes a uniform distribution of emis- 
sion over 2° x 2° on the sky, the CXB model for 4 square degree 
field in units of photons crrr^-'keV- 1 is given by 



— =9.412xl0~ 3 exp . (1) 

dE UkeV/ \ 40keV/ 



The 20-50 keV energy flux of the CXB model (Eq. Q}, 
5.8 x 10 8 erg s 1 cm 2 sr 1 , agrees with the recen t report from 
the Beppo-SAX observations dFrontera'e t al. 2007) within about 
8%. This level of uncertainty in the CXB model is smaller than 
the systematic error of the current NXB model. However, ac- 
cording to the calibration report based on the Crab observations, 
the PIN spectral normali zation is systematically higher than the 
previous results given in iToor & Sewarl ( fl97l by 13%. Thus, 
we increase CXB given in Eq.[T]by a factor of 1.13 in the simu- 
lation. 

The resulting simulated CXB spectrum is shown in Fig. [7] 
To calculate the CXB spectrum the integration time of 100 Ms 
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Fig. 7. Observed PIN spectrum minus the nominal NXB model 
(Data-NXB; black), the simulated CXB model (red), and the 
point source contribution estimated from the XMM-Newton ob- 
servations of the same field (green). 
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Fig. 8. Observed PIN spectrum with the background model sub- 
tracted (Data-NXB-CXB; the crosses) in the 12-60 keV band. 
The error-bars show only the statistical errors. The histogram 
shows the best-fitting APEC thermal plasma model. The bottom 
panel shows the residuals in units of the number of standard de- 
viations. 



was assumed; thus, the statistical error of the simulated spectrum 
is negligibly small. The CXB count rates in the 12^40 keV and 
40-60 keV bands are given by 2.9 x 10~ 2 and 6.1 x 10~ 4 c s -1 , 
respectively, which amount to about 6.1% and 1 .2% of the NXB, 
respectively. 

Kushino et al. (2002) have found that, based upon the ASCA 
observations, the CXB fluctuation in the 2-10 keV band is ~ 6%, 
which can be attributed to the Poisson noise of the source count 
within 0.5 degree 2 . Moreover, by scaling the HEAO-1 result of 
2.8% with the equation ocxb <* Q~ ()5 S ()25 (Q and S are the 
effective solid angle of the observation and upper cut-off flux 
of detectable discrete sources in the field of view, respectively; 
ICondodfl974l) . iKawano et al.l d2008l) calculated the lcr fluctua- 
tion of CXB to be 9.2% for the ~ 100 ks HXD/PIN observations. 
Therefore, the CXB fluctuation over the 30' scale is small com- 
pared to the NXB systematic error. 

4.3. Hard X-ray point sources 

Next, we estimate the hard X-ray flux of point sources inside the 
field of view of HXD/PIN from the 30 ks XMM-Newton/PN data 
on the same field (Observation ID: 0112960101). The data re- 
duction was carried out in the standard manner with SAS version 
7.1.0 and the periods of background flares were removed with 
the threshold value of 0.22 c s _1 in the 12-14 keV band. 

With the edetect_chain program, thirty point sources have 
been detected in the PN image. The photon index for each 
source was estimated from the hardness ratio, defined as a ra- 
tio of the source count rates in the 0.5-1.5 keV and 1.5-5 keV 
band. The backgroun d was subtracted using the blank-sky data 
dCarter & Read 2007). The galactic absorption was assumed in 
the calculation. The photon index of F > 2 was obtained for al- 
most all the sources that have significant emission in the hard 
band. We then simulate the PIN spectrum expected for each 
source using a power-law spectrum with V - 2, where we in- 
clude the PIN angular response function using hxdarfgen. 

Fig.Qshows the spectrum of the sum of the detected sources. 
We find that the contribution from the point sources is totally 
negligible. 



4.4. PIN source spectrum 

In Fig. [H the observed PIN spectrum with NXB and CXB sub- 
tracted is shown. In the 12-40 keV band we detect the hard X-ray 
emission at the 9cr level. Taking into account the NXB system- 
atic error, however, significance of the detection becomes much 
lower: the estimated countrate is 0.0187±0.0021 (±0.0142)c s _1 
in the 12^40 keV, where the first and second errors are the 
lcr statistical and lcr systematic uncertainties, respectively. At 
higher energies, 40-60 keV, the count rate is 7 ± 6 (+15) x 
10~ 4 cs~' (lcr); thus, we do not expect significant hard X-ray 
emission in the 40-60 keV band. 

Next, we fit the single-temperature APEC model to the 
background-subtracted PIN spectrum in the 12-60 keV band. 
(The XIS data is not used here. For a joint analysis of XIS and 
HXD/PIN, see §0) The metal abundance was fixed at the best- 
fitting value obtained from the XIS analysis, 0.33 solar. The sys- 
tematic error of the spectral parameters due to the systematic 
uncertainty in the NXB model was estimated by changing the 
normalization factor of the NXB model by +3% and repeating 
the fitting procedure. 

The results are shown in Fig. |8]and Table|6] We find that the 
temperature and normalization of the single-temperature APEC 
model are kT = 20.2+ 21 4 (+£?,) keV and 1 .02+°^ (+° ^) x 10~ 2 , 
respectively. These parameters are consistent with those ob- 
tained from the XIS data. Therefore, we conclude that the hard 
X-ray emission observed with PIN does indeed come from 
RXJ1347.5-1145. 

From the observed PIN count rate and the best-fitting single- 
temperature APEC model given in Table [2] we estimate the en- 
ergy flux to be 3.8 ± 0.5 (±3.3) x 10~ 12 erg s^'cirr 2 in the 12- 
60 keV (lcr). On the other hand, the previous Beppo-SAX obser- 
vation of RX J1347.5— 1 145 has reported detection of the hard 
X-ray emission at the ~ 1 . 5cr significance: (4 .0 ±2.6) x 10~ 2 c s _1 



in the 13-60 keV band dEttori et al.1 [2001b. and (6.8 ± 4.4) x 
10~ 2 c S -1 in the 20-80 keV band (Nevalainen et al. 2004). With 



the 12.9 keV Raymond-Smith model and the response of the 
PDS instrument on board Beppo-SAX, the pimms program gives 
the 12-60 keV flux as (4.5 ± 2.9) x 10~ 12 erg s cirr 1 . (All errors 
quoted here are lcr.) Thus, the PIN result agrees with that of 
PDS to within the lcr errors. 
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Table 6. Single-temperature APEC model parameters fit to the PIN data in the 12-60 keV band. The systematic error in the NXB 
model is included in the error estimate. 



Sensor 


kT [keV] 


Z [solar] 


K 


X z /dof 


NXB a 


PIN 


20.2(12.2-41.6) 


0.33(Fix) 


1. 02(0.59-1. 97)x 10-" 


6.7/8 


1.00 


PIN 


11.5 


0.33(Fix) 


0.76 xl0~ 2 


26.2/8 


1.03 


PIN 


27.9 


0.33(Fix) 


1.28xl0- 2 


6.9/8 


0.97 



a The NXB normalization factor. 



5. XIS+HXD joint analysis 

In §[3]the Suzaku XIS data indicated the presence of a significant 
amount of very hot gas. In this section, we study the physical 
properties of this hot component, such as the amount, tempera- 
ture, and metal abundance, by performing a joint analysis of the 
Suzaku data and the spatially-resolved Chandra spectra. 

In § 15. li the properties of the ambient gas are derived from 
the Chandra data. In § l5.2l the properties of the hot gas compo- 
nent in the SE quadrant are derived from the Suzaku broad-band 
spectra combined with the Chandra data by means of a multi- 
temperature model. In § 15.31 we argue that the hot component 
is explained better by thermal gas than by non-thermal gas, and 
derive an upper limit on the possible non-thermal emission. 

5. 1 . Modeling the cluster average component with Chandra 

The deep Chandra ACIS-I data on RX J1347.5-1145 (Obs 
ID:3592, Date:2003-09-03) was analyzed with CIA0 version 3.4 
and CALDB version 3.3.0. The net exposure time after removing 
the periods of high background rates is 56.1 ks. The backgrounds 
were estimated from the same detector regions of the blank-sky 
data, whose normalization factors were determined with the ra- 
tios of the 10-12 keV count rates. 

Fig. [9] shows the Chandra image and the azimuthally- 
averaged surface brightness profiles for three different position 
angles (0° - 360°, NW:-90° - 180°, and SE:180° - 270°). The 
cluster emission extends out to r = 5' ~ 1.7 Mpc and the sig- 
nificant excess emission is visible in 10" < r < 60" for the SE 
quadrant. 

To derive the average temperature profile within 5', we mea- 
sure the spectra in the following six radial bins in the NW re- 
gion, i.e., the region outside of the SE quadrant: - 4", 4" - 8", 
8" -24", 24" - 56", 56" - 1 20", and 1 20" - 300" . The radial bins 
were chosen such that the temperature profile is measured with 
sufficient spatial resolution and the statistical accuracy ~ 20%. 

Fig. [10] shows the projected temperature profile, i.e., the 
temperature profile determined from the simple APEC model 
fit to each radial bin in the NW region. Our estimate of the 
projected temperature profile derived from the Chandra ACIS- 
I data ag rees with the pre vious Chandra ACIS-S results re- 
ported by Alle n et al. 1 (120021) to within the statistical error, partic- 
ularly in the central 100"region. We have also confirmed based 
on the ACIS-I spectral simulations assuming the same photon 
statistics as the present observation that we can securely deter- 
mine the temperature of each bin including the high-temperature 
(~ 18 keV) regions between 24"and 120". 

Since each spectrum is a superposition of the cluster emis- 
sion from any points along the line of sight, we correct this ef- 
fect, or "deproject," by using the PROJCT model in the XSPEC 
software and fitting the six radial bins simultaneously under the 
assumption that the temperatur e distribution is spherically sym- 
metric. As noted in lAllen etal.ld2002l) . the metal abundance gra- 
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Fig. 10. (a) Projected (crosses) and deprojected (circles) temper- 
ature profiles measured from the Chandra data in the NW region, 
(b) Projected metal abundance profile (crosses). The horizontal 
dashed and dotted lines show the best-fitting constant abundance 
and the corresponding 90% confidence interval obtained from 
the deprojection analysis, respectively. We shall assume that the 
metal abundance is constant over all radial bins throughout this 
paper. 



dient is marginally seen in the Chandra data; however, the large 
uncertainty in the metal abundance measurement does not war- 
rant our treating it as a free function. We therefore assume that 
the metal abundance distribution is uniform and ignore the abun- 
dance gradient. 

The resulting deprojected temperature profile in the NW re- 
gion is shown in Fig.[10]and the best-fitting parameters are given 
in Table [7] Henceforth, we refer to the best-fitting model pre- 
sented in Table|7]as the "6APEC model." 
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Fig. 9. (a) Chandra ACIS-I image in the 0.5-7 keV band. The solid lines show the NW and SE regions. The radius of the circle is 
5'. (b) The surface brightness profiles of the position angle 0° - 360° (black), NW:-90° - 180° (red), and SE:180° - 270° (green). 
The black line shows the best-fitting double yS-model for the position angle of 0° - 360°. 

Table 7. Parameters of the "6APEC model," estimates of the deprojected temperature profile from the Chandra ACIS-I data at the 
six annular bins in the NW region, which excludes the SE quadrant. 



Radius ["] 


kT [keV] 


Z [solar]" 


K b 




x'/dof 


0-4 


5.74 (5.01 -6.64) 


0.47 (0.39-0.55) 


0.97(0.92 - 


1.02) x 10 -3 


364.2/290 


4-8 


8.06 (7.10-9.20) 




1.83(1.77- 


1.89) x 10~ 3 




8-24 


14.12(12.61 - 16.08) 




4.02(3.93 - 


4.10) x 10" 3 




24-56 


17.49 (14.95-21.26) 




2.86(2.79 - 


2.95) x 10~ 3 




56-120 


18.47 (14.49-23.18) 




2.29(2.22 - 


2.37) x 10" 3 




120-300 


8.27 (6.15-13.16) 




0.91(0.86- 


0.97) x 10~ 3 





a The metal abundance is assumed to be common for all radial bins. 

b K represents the normalization factor for each spherical shell, i.e., the position angle of 0° - 360°, of the fitting. 



5.2. Temperature measurement of the SE clump 

The 6APEC model in § 15.11 describes the average temperature 
structure of ambient gas, measured from the NW region. We sub- 
tract the flux accounted by this model from the Suzaku data, and 
study the nature of the excess emission in the SE quadrant. 

We add another APEC model to describe the hot component 
in the SE quadrant. We use the best-fitting parameters for kT 
and K given in Table [TJfor the ambient gas (6 APEC model), and 
use the Suzaku data, both the XIS and PIN spectra, to constrain 
the temperature, kT ex , and the normalization, K ex , of the excess 
component in the SE quadrant. As for the metal abundance, the 
abundance is assumed to be equal to all of the six radial bins of 
the ambient gas, as well as to the excess component, i.e., a single 
abundance is applied to the entire cluster. We do this because the 
metal abundance is not constrained by the Chandra data very 
well. 

We have found that the Chandra data tend to give system- 
atically higher values for kT and K than the XIS data do. In 
order to correct for the difference between their instrumental cal- 
ibrations, we multiply the 6APEC model by the PLABS model, 
M(E) = E~ a , with a = 0.02 (E is the energy of photons), and 
use the relative normalization factor of 0.93. See AppendixlAlfor 
more details on the Suzaku-Chandra cross-calibration. 

Fig. [Hand Table[8]show the results of the 6APECxPLABS 
+ APEC model fit to the XIS (0.5-10 keV) and PIN (12-60 keV) 



data. To show the degree of systematic errors due to the uncer- 
tainty in the Suzaku-Chandra cross-calibration, we also show 
the parameters for a = and 0.03 in Table [8] From these re- 
sults we conclude that the temperature of the excess emission 
in the SE quadrant is kT ex = 25.3+^ (+^) keV, where the 
first error is statistical and the latter is systematic (both 90% 
C.L.). In Table [8] we also show the fitted parameters when only 
XIS was used. Thanks to the high photon statistics of XIS, the 
temperature of the SE excess component can be constrained 
as kT™ s = 24.8^-g Cf s 5 ) keV. However, the spectral cut-off at 
25/(1 + z) ~ 17 keV actually falls into the PIN band, which 
suggests the importance of the broad-band spectroscopy in the 
present study. As for the other quantities, we find from the 
XIS+PIN joint analysis that the metal abundance is 0.38(0.35- 
0.41) solar, and the absorption-corrected X-ray luminosity of the 
excess component is Lx.ex = 1.4xl0 45 ergs -1 in the 0.5-10 keV, 
and 1.2 x 10 45 erg s 1 in the 12-60 keV. 

Since the determination of kT ex from the X-ray spectroscopy 
is the central result of this paper, we estimate the other potential 
source of systematic errors in our analysis. Specifically, we ex- 
amine: (i) the assumption of constant metal abundance over all 
radii, (ii) some arbitrariness in the choice of the SE region, and 
(iii) effects of the uncertainty in the 6APEC model parameters. 

(i) What if the metal abundance of the ambient gas and that 
of the excess component are different? What if the abun- 
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Table 8. Spectral parameters for the SE excess component obtained from the 6APECxPLABS+APEC model fit to the XIS+PIN 
data and the XIS data. 







APEC 




* 2 /dof 




Sensor 


kT, x [keV] 


Z ex [solar]" 




a b 


XIS0-3+PIN 


25.3(20.8-31.4) 


0.38(0.35-0.41) 


2.14(2.05-2.24)xl0- j 


1311.1/1219 


0.02 


XIS0-3+PIN 


15.8 


0.36 


1. 98x1 3 


1313.1/1219 


0.00 


XIS0-3+PIN 


32.2 


0.39 


2.26xl0- 3 


1310.2/1219 


0.03 


XISO-3 


24.8(20.2-30.9) 


0.38(0.35-0.41) 


2.00(1.91-2.09)xlO- J 


1303.3/1209 


0.02 


XISO-3 


15.5 


0.36 


1.84X10 3 


1303.7/1209 


0.00 


XISO-3 


31.9 


0.39 


2.11X10- 3 


1303.0/1209 


0.03 



a The metal abundance is assumed to be the same for all the APEC components, including the ambient gas as well as the excess component. 
b The power-law index of the PLABS model, coming from the Suzaku-Chandra cross-calibration (see text). 



dance gradient does exist in the ambient gas? To explore 
these possibilities we treat the metal abundance of the ex- 
cess component, Z ex , as well as that of the ambient gas in 
two larger radial bins, 0" < r < 8" and 8" < r < 300", 
as free parameters, and fit them simultaneously. We find 
kT ex = 26.3 keV, Z ex < 1.8 solar, and K ex = 2.16 x 10~ 3 
(X 2 /dof = 1306/1217), which are consistent with the values 
given in Table [8] As for the metal abundance of the ambient 
gas we find Z ~ 0.64 solar and 0.22 solar in 0" < r < 8" and 
8" < r < 300", respectively. 

(ii) From the Chandra image one may observe that the excess 
emission appears to extend over slightly more than the SE 
quadrant. How does this affect the determination of kT ex l We 
repeat our analysis with the SE and NW regions re-defined 
as (180° - 315°) and (-45° - 180°), respectively. We find 
kT ex = 20.0keV,# ex = 2.87 xl0~ 3 , andZ ex = Z = 0.37 solar 
with ;f 2 /dof = 1315/1219. Therefore, the exact choice of the 
SE region does not change the best-fitting parameters very 
much. 

(iii) The statistical errors in the 6APEC model parameters had 
not been propagated through the spectral fitting of the SE 
quadrant, except for the cross-calibration error in the spec- 
tral slope, which was included by varying a from to 0.03. 
To address this, we carry out simultaneous APEC fits to 
the Chandra spectra in 6 radial bins (without deprojection) 
as well as to the Suzaku XIS and PIN spectra, for which 
we use the 6 APEC models multiplied by the a = 0.02 
PLABS model plus the additional APEC model. We find 
kT ex > 26 keV, Z ex = 0.38(0.35 - 0.41) solar, and K ex = 
2.16(1.96 - 2.41) x 10~ 3 (x 2 /dof = 1676/1510). Although 
the temperature is not well constrained in this case, it equally 
shows a high value. 

From these studies we conclude that our estimate, kT ex = 
25.3^-1 (+6-9) keV, is robust. 

Finally, as a consistency check we analyze the Chandra 
spectrum of the SE quadrant (r < 5') on its own, without 
the Suzaku data, in the same manner. The best-fitting parame- 
ters for the excess component are kT^ x andm = 31.1*^H keV, 
z chandra _ 0.49+J}^ solar, and Kg hmdra = 1.47+°^ x 10~ 3 
(X 2 /dof = 320/332). Therefore, both the Chandra-alone analy- 
sis and the joint Chandra+Suzaku analysis show that the excess 
emission has the temperature in excess of 20 keV0 



5 As for the normalization factor, K^"" dm is found to be smaller than 
that from the joint analysis by about 35%. This may be explained by the 
difference in their integration regions. On the other hand, this level of 
difference actually amounts only to a few% of the total cluster emis- 



Compared to the results from the Chandra data alone, the 
joint Chandra+Suzaku broad-band data analysis yields much 
more accurate determination of kT ex . The Suzaku's unprece- 
dented sensitivity over the wide X-ray band makes it possible 
to determine, for the first time, the temperature of such hot gas 
in the ICM solely from the X-ray spectroscopy without help of 
the SZ effect. 

5.3. Constraint on the non-thermal emission 

So far, we have been assuming that the excess hard X-ray emis- 
sion from the SE quadrant is thermal; however, could it actually 
be non-thermal? If so, our derived temperature must be recon- 
sidered. 

To address this question, the XIS and PIN data is re-analyzed 
in light of non-thermal emission. While the 6APEC thermal 
plasma model is again used for the ambient gas, the APEC 
model is replaced with a power-law (PL) spectrum for describ- 
ing the excess component. We then fit the new model, "6APECx 
PLABS +PL model," to the observed XIS and PIN data. The 
amplitude (normalization), a power-law index, F, of the non- 
thermal component, and the metal abundance of the ambient gas 
are treated as free parameters. 

Fig. [12] shows the results, and Table [9] lists the best-fitting 
parameters. We find that the power-law index is given by F ~ 
1.5, or more precisely F = 1.45*2'^ Cqo^)- The reduced^ 2 of 
this model, x 2 = 1317.1/1219 (see Tabled, is slightly larger 
than that for the thermal model, x 2 = 1 3 1 5 . 1 / 1 2 1 9 (see Table©, 
and negative residuals are seen in consecutive bins above 22 keV. 

Let us examine the non-thermal model more closely. In the 
XIS band below 10 keV, the non-thermal model with F = 1 .5 and 
the thermal APEC model with 25 keV are hardly distinguishable. 
However, the observed PIN spectrum is much softer than T = 
1.5, with the effective photon index being r e fr ~ 3. When we 
fit the 6APECX PLABS+PL model to the PIN data alone, the 
90% lower bound on the photon index is found as F > 1.83, 
which is significantly above the 90% upper bound on F from 
the XIS data, F < 1.5. We have further checked that fixing the 
power-law index to F = 2.0, as predicted f rom the non-thermal 
bremsstrahlung in the strong shock limit dSarazin & Kem pner 
2000), yields a poor fit to the XIS and PIN data; the reduced^ 2 is 
1708 for 1220 degrees of freedom. Therefore, it seems difficult 



sion, and is comparable to the calibration errors of the instruments. 
Therefore, we simply assign 30% systematic error to A* cx . 

6 The systematic errors due to both the NXB model and the Suzaku- 
Chandra cross-calibration are included. 
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Table 9. Spectral parameters for the SE excess component obtained from the 6APECxPLABS+non-thermal power-law (PL) model 
fit to the XIS+PIN data. 







Power-law 


6APEC 


* 2 /dof 




Sensor 


r 


Normalization 


Z [solar] 


a 


XIS0-3+PIN 


1.45(1.41- 


-1.48) 2.58(2.30-2.52)xlfr 4 


0.42(0.39-0.45) 


1317.1/1219 


0.02 


XIS0-3+PIN 


1.54 


2.57X10- 4 


0.40 


1312.5/1219 


0.00 


XIS0-3+PIN 


1.41 


2.58xl(T 4 


0.42 


1321./1219 


0.03 



Table 10. Limits on the amplitude (normalization) of a non-thermal power-law (PL) component with F = 1 .5 as the explanation of 
the excess hard X-ray emission in the SE quadrant. 



Sensor 


Model 3 


Normalization" 




^r/dof 


Flux(12-60keV) d 






NXBxl.00 


NXBx0.97 


NXBxl.00 


NXBxO.97 


[ergs -1 cm" 2 ] 


PIN 


APEC+PL 


~ 


5.1 x 10 -* 


7.2/7 


4.6/7 


< 2.11 x 10 -11 


XIS0-3+PIN 


APEC+PL 


< 1.6 x 10" 4 


2.0 x 10~ 4 


1310.2/1218 


1344.3/1218 


< 8.2 x 10~ 12 


PIN 


PL 


1.8(0.7 - 2.9) x 10-" 


6.3 x 10~ 4 


8.8/9 


5.1/9 


<2.15x 10-" 


XIS0-3+PIN 


PL 


2.48(2.39 - 2.59) x 10~ 4 


2.52 x 10~ 4 


1322.8/1220 


1351.9/1220 


<3.8x 10~ 12 



a The spectral model for the SE excess emission. 

b The amplitude (normalization) of the power-law component after removing the NXB models with rescaling factors of 1 and 0.97. The error 

and limits are the lcr bound. 
c The values of x 1 and the number of degrees of freedom. 

d The 3cr upper limit on the energy flux of the power-law component in the 12-60 keV band. 



to explain the SE excess component with the non-thermal PL 
model. 

Given that the excess component is explained better by the 
thermal emission than by the non-thermal emission, we derive 
an upper limit on the non-thermal emission as the explanation 
for the SE component. We fix the metal abundance and the PL 
index at 0.38 solar and r = 1 .5, respectively. The 6APEC model 
modified by the a - 0.02 PLABS model is again used for de- 
scribing the ambient gas. 

The derived upper limits are given in Table [10] Taking into 
account the NXB systematic error of the PIN detector, the 3<x 
upper limit on the non-thermal flux in the 12-60 keV, Fhxr, 
is estimated as Fhxr < 2.1 x 10~ n ergs^'cm -2 from the PIN 
analysis, and F H xr < 8 x 10~ 12 ergs _1 cnr 2 from the XIS+PIN 
simultaneous fit. Here, lcr error is calculated by adding the lcr 
statistical error and the lcr systematic error in quadrature. 

6. Discussion 

Using the deep broad-band observations of RX J1347.5— 1 145 
with Suzaku, the temperature structure has been determined with 
a particular focus on the extremely hot gas. We have found 
that the multi-temperature thermal emission model including the 
very hot (kT ex ~ 25 keV) component fits the wide-band spectra 
well. We have also placed an upper limit on the non-thermal hard 
X-ray emission. 

In § 16.11 our results are compared with the previous X-ray 
and radio observations, and examine the physical properties of 
the SE clump further. In § 16.21 our results are compared with 
the radio mini halo in this cluster that has been discovered from 
the 1.4 GHz radio observations to estimate the magnetic field 
strength in the ICM. 

6.1. Properties of the extremely hot gas 

From the analysis of the Suzaku spectra of RX J1347.5-1 145, 
the temperature of the SE excess emission is obtained to be 
kT ex = 25.3+^ C 6 9 9 5 ) keV (90% C.L.; the first error being sta- 



tistical and the second systematic). This is in an excellen t agree- 
ment with the previous measurement by Kitava ma et al.l {2004), 
28.5 + 7.3 keV (68%; statistical only), from a joint analysis of 
the SZ d ata dKomatsu et alj 1 999l l200lb and the Chandra ACIS- 
S3 data (Alle n et al. 2002 ). Here, we compare our X-ray-only 
results with their SZ+X-ray results. 

We estimate the gas density and gas mass of the SE excess 
component from the Suzaku and Chandra data. For simplicity, 
we assume that the extremely hot gas is uniformly distributed 
within a sphere of a radius R - 25" ~ 144 kpc (because the ex- 
cess is present in 10" < r < 60"; Fig- |9j>). Using the measured 
spectral normalization, K ex , which depends on n e nuV (where 
V = 4jtR 3 /3 is the volume of the gas and «h = 0.86n e ), the 
electron density, n ex , and the gas mass, M gas ex , are obtained as: 



n e , ex = (1.6 ± 0.2) x 10- 2 ( — V " cm" 3 , (2) 

\3.6x 10 71 cm 3 / 

M gas , ex = (5.6±0.8)xl0 12 ( ^» ,) M Q . (3) 

V3.6 x 10" cm J / 

Here, we have propagated the systematic error in K ex through 
the final results (§ |5_^2j> . These r e sults may be compared with 
those obtained in Kitav ama et all d2004l) : n ex = (1.49 + 0.59) x 
10~ 2 cirr 3 and M gaSjex ~ 2 X 10 12 M©. While the inferred num- 
ber densities are in an excellent agreement, the gas mass from 
our analysis is more th an a factor of two greater than that from 
iKitavama et alj d2004l); however, the uncertainty in the gas mass 
estimate from Kitavama et alj d2004l) is large enough for them to 
be consistent. (The uncertainty is large because they did not as- 
sume spherical geometry for the excess component, and the line- 
of-sight extension of the SE component was poorly constrained.) 

How does the excess component compare with the rest of the 
cluster? The average (ambient) temperature and gas density for 
10" < r < 60", excluding the SE quadrant, are kT ~ 15 keV 
from the C handra spectrum and n e ~ 6.6 x 10~ 3 cirr 3 from 
Eqs. 3-5 of IKitavama et al. (2004), respectively. Thus, the SE 
clump exhibits the temperature and the density that are higher 
than the ambient gas in the same radial bins by factors of 1.6 
and 2.4, respectively. This means that the excess hot component 



N. Ota et al.: Suzaku broad-band spectroscopy of RX J1347.5— 1 145 



13 



> 




Energy [keV] 



~ (b) 




Energy [keV] 

Fig. 11. (a) 6APECxPLABS(a = 0.02) + APEC model, fit 
simultaneously to the observed XIS (0.5-10 keV) and PIN 
(12-60 keV) spectra. The spectra taken by four sensors, XIS- 
(black), 1 (red), 2 (green), and 3 (blue), are shown sepa- 
rately. The black and red lines show the ambient component 
(6APECxPLABS) and the excess hot component described by 
an additional APEC model, respectively, convolved with the 
telescope and the detector response functions. In the bottom 
panel the residuals of the fit in units of the number of standard 
deviations are shown. We have used the nominal NXB model for 
the background of the PIN data (see § I4.U . (b) same as (a), but 
for the unfolded XIS-0 and PIN spectra. 



> 




Energy [keV] 



- (b) 



> 1 




Energy [keV] 

Fig. 12. (a) 6APECxPLABS(a = 0.02) (thermal model for the 
ambient gas) + PL model (non-thermal model for the hot gas), 
fit simultaneously to the observed XIS (0.5-10 keV) and PIN 
(12-60 keV) spectra. The spectra taken by four sensors, XIS- 
(black), 1 (red), 2 (green), and 3 (blue), are shown sepa- 
rately. The black and red lines show the ambient component 
(6APECxPLABS) and the excess hot component described by 
a non-thermal power-law model, respectively, convolved with 
the telescope and the detector response functions. In the bottom 
panel the residuals of the fit in units of the number of standard 
deviations are shown. We have used the nominal NXB model for 
the background of the PIN data (see § 14. Il l, (b) same as (a), but 
for the unfolded XIS-0 and PIN spectra. 



is over-pre ssured, and such a region is expected to be short-lived, 
~ 0.5 Gyr dTakizawall 19991) . 

As already discussed in lKitavama et al. (2004), the gas prop- 
erties can be explained by a fairly recent (within the last 0.5 Gyr 
or so), bullet-like high velocity (Av ~ 4500 kms -1 ) collision of 
two massive (5 x 1O 14 M0) clusters. 

Moreover, the heat energy of the SE clump is estimated to 
be Eth,ex = (3/2)AT„(7i wat + « H ,ex)V ex ~ 6 x 10 62 erg. This huge 
amount of energy cannot be easily produced by a central source 
in the cluster: one would need an AGN with 10 46 erg s for at 
least 1 Gyr and put the energy into the SE region without any 
radiative loss. On the other hand, cluster mergers, which are the 
most energetic events in the Universe after the Big Bang, will 
most naturally explain this high energy phenomenon. Therefore 
our results do support the merger scenario, solely from the X-ray 
spectroscopic data without help of the SZ data. 



6.2. Estimation of the magnetic field 

What physics can we learn from our upper limit on the non- 
thermal hard X-ray emission from the Suzaku HXD/PIN? The 
non-thermal hard X-ray emission is produced via the inverse 
Compton (IC) scattering of relativistic electrons off the Cosmic 
Microwave Background (CMB) photons, and the same popula- 
tion of electrons also p roduce the synchrotron radia tion. From 
the exact derivations by Blumenthal & Gould ( 1970), equations 
for the synchrotron emission at the frequency vs yn and the IC 
emission at vie are: 



dW s 



4nN e 3 B ( P +1)/2 I 3e 



dv Syn dt 

dW lc 
dvicdt 



4mn„c 



-<7>-l)/2 



a{p>l%- l)l \ (4) 



^/z-C +3 >/^o(^r CMB )^ 5 )/ 2 F (/ ,)v^- 1)/2 . (5) 
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Here No and p are the normalization and the power-law index of 
the electron distribution, N(y) - Noy~ p (y is the Lorentz factor 
of the electron), h is the Planck constant, 7cmb is the CMB tem- 
perature and 7cmb = 2.73(1 + z) K. The functions a(p) and F{p) 
are defined as follows: 



a(p) 
F(p) 



2 (p-D/2V3r(^i)r(^)r(^ 



8^/2 (p+ i )r (£±Z) 



2P +Z {p l +4p+ 11) r(P + 5 \ (P + 5 



(p + 3) 2 (p + 5)(p + l) 



(6) 
(7) 



Eqs.H] and |5] give the total emitted powers per volume per fre- 
quency in the rest frame. Given that both powers are diminished 
by the same dimming factor, the ratio of observed flux densi- 
ties of the synchrotron and IC emission, Ssyn/Sic, is equal to 
(dWsyn/dvs yn dt)/(dWic/dvicdt), and the strength of the mag- 
netic field B can be directly estimated. 

We find S\c < 0.11 p]y from our limit on the non-thermal 
hard X-ray emission for electrons with the IC emission energy 
(frequency) of 12 k eV (v IC = 2.9x 10 18 (1 +z) ' Hz). As forS S y„, 
iGitti e t al. (2007a) discovered extended (~ 500 kpc) radio emis- 
sion in RX J1347.5-1145 based on their low (18") resolution 
VLA observations and higher (2") resolution data in the VLA 
archive. They subtracted contribution of discrete radio sources 
and estimated the total flux density of the diffuse radio emission 
to be Ssyn = 25 mJy at vs yn = 1.4(1 + z) _1 GHz. Combining 
these numbers with Eqs. © and (0, a lower bound on the mag- 
netic field strength in the ICM is obtained as B > 0.007 pG for 
p = 2Y — 1 = 2. The estimation is sensitive to the assumption of 
p: for example, in the case of p — 3, B > 0.072 pG. 

This limit, though weak, is consistent with typical values 
found in other clusters, B ~ 0.1 — 1 nG, based on the RXTE 
and Beppo-SAX observations (e.g., Rep haeli et alj2008T) . Our re- 
sult is also comparable to the recent Suzak u report on a mergin g 
cluster A3376 at z = 0.046, B > 0.03 pG dKawano et alj |2008). 
However, the previous measurements came mostly from nearby 
(z < 0.1) clusters as well as from some medium-r edshift clusters 
such as A2163 at z = 0.20 dReohaeli et al.ll2006l) and the Bullet 
cluster at z — 0.296 (Petrosia n et al.ll2006l) . Our work provides a 
constraint on the cluster magnetic field strength at a higher red- 
shift, z = 0.45 1 . 



7. Summary 

We have reported on the results from the analysis of the 
Suzaku wide-band (0.5-10 keV with XIS and 12-60 keV with 
HXD/PIN) X-ray spectroscopic observations of the most lumi- 
nous X-ray cluster, RX J1347.5-1 145, at z = 0.451, in order to 
investigate the temperature structure of the ICM, signatures of a 
recent violent merger, as well as signatures of the non-thermal 
emission. 

This cluster is known to contain a very hot gas clump that 
produces the excess emission on top of the a mbient gas that 
follows more-or-less the conventional B mod el dKomatsu et alj 
l200UlAllen et al.ll2002t iKitavama et al.ll2004l) . The temperature 
of this gas clump has not been measured accurately from the 
previous X-ray spectroscopic observations, as their sensitivities 
degrade significantly beyond 10 keV. 

The re-analysis of the Chandra data confirms the previous 
work, and yields a poor limit on the temperature of the excess 
hot component, kT^ wndra = 31.1^ keV (90% C.L.; statisti- 
cal). When the hard X-ray data from the Suzaku XIS spectra in 
the 0.5-10 keV region and the HXD/PIN data in the 12-60 keV 



are combined with the Chandra data, we finally obtain a good 
measurement of the temperature, 25.3;^ C.f 9 5 ) keV (90% C.L.; 
statistical and systematic), which is in an excellent agreement 
with that derived from the previous joi nt analysis of the S Z ef- 
fect and X-ray imaging observations (IKitavama et al.1120041) . We 
stress that this is the first time that the X-ray spectroscopic obser- 
vations alone give a good handle on such a high temperature gas 
component, which is made possible by Suzaku s unprecedented 
sensitivity over the wide X-ray band. 

We have found that the broad-band X-ray spectrum is ex- 
plained better by the thermal plasma model than by the non- 
thermal power-law model. Thus, the present result confirms the 
presence of the hottest thermal gas in the cluster. The most likely 
explanation for this phenomena is a recent violent merger with 
the collision velocit y of ~ 4500 km/s, similar to the one found 
in the Bullet clusterdMilosavlievic et alj2007 : Spring eT& Farrarl 
l2007tlMastropietro & Burkertl2007tNusserfl2008l) . 

The upper bound on the non-thermal flux in the 12-60 keV 
band, Fhxr < 8 x 10~ 12 ergs^'crrT 2 , yields, when combined 
with a recent disco very of the radio mini halo in this cluster by 
IGitti et aLl (l2007al) . a lower bound on the ICM magnetic field 
strength, B > 0.007 pG for the electron index of 2. These con- 
straints provide valuable information on the non-thermal nature 
of the cluster gas and the particle acceleration in a high-redshift 
universe. However, the accuracy of the present measurement of 
the hard X-ray emission is limited by the sensitivity of the in- 
strument. Suppose that the magnetic field in RX J1347.5— 1 145 
is comparable to that of the other clusters, B ~ 0.1 pG, a 50-fold 
higher sensitivity is required to detect non-thermal hard X-ray 
emission from this cluster. We expect th at the future X-ray mis- 
sions (e.g., NeXT: \Takahashi et al.l2008l) will determine the non- 
thermal X-ray flux accurately and draw a more complete picture. 

Before the advent of the Suzaku satellite, it was difficult to 
measure the high-temperature (» 10 keV) gas in the X-ray band 
reliably, due to the limited sensitivity. The present paper demon- 
strates the power of the broad-band X-ray spectroscopy in the 
study of gas heating due to the cluster merger, which is nicely 
complementary to the observations of the SZ effect. In the near 
future we expect the number of samples of merging clusters ob- 
served with Suzaku to increase rapidly, which would shed new 
light on the physics of violent mergers. 
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Appendix A: The Suzaku-Chandra cross calibration 

Because the temperature measurement of very hot gas based on 
the spectral fit is sensitive to a subtle change of the spectral slope, 
a precise calibration of the effective area is critical. In the analy- 
sis of XIS+HXD spectra (§|5), we have modeled the cluster aver- 
age emission with Chandra. Thus the cross-calibration between 
Suzaku/XIS and Chandra/ ACIS-I is mandatory. 

If the global cluster spectrum is accumulated from the r < 5' 
circular region with the Chandra/AClS-l data, the APEC model 
fitting gave kT = 14.4(13.6 - 15.3) keV, Z = 0.45(0.38 - 0.52), 
K = 1.47(1.45 - 1.49) x 10~ 2 Cr 2 /dof = 552.8/426), where 
Z = 0.451 and jV H = 4.85 x 10 20 ctrT 2 are adopted. In compar- 
ison with the best-fit XIS parameters (Table 0, the temperature 
and the normalization are higher by about 1.5 keV and 7% for 
Chandra. We discuss below the possible reasons in the light of 
the Suzaku and Chandra calibration uncertainties. 

1 . Suzaku XIS 

Since the spatial extent of X-ray emission from RX J1347.5- 
1145 is small compared to the PSF of Suzaku, it can be 
treated as a point source (we have confirmed that the spec- 
tral fit with an auxiliary file generated for a point source 
at the HXD-nominal position yields statistically consistent 
spectral parameters with those listed in Table [2j. According 
to the latest calibration of XIS and PIN, the best-fit power- 
law parameters describing the Crab spectra are the photon 
index F = 2.073 ± 0.006 and the normalization 9.21 ± 
0.1 photons cm _2 s _I keV _1 , which are by about 1.5% and 5% 
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Fig. A.l. Suzaku/XIS + Chandra/ ACIS -I spectra of RX J 1347.5- 
1145 (r < 5') simultaneously fitted with the APECxPLABS 
model 



lower than the 'standard' val ues (T = 2.10, the normalization 
= 9.7: lToor&Sewardlll974 . 

To simply simulate the effect of T, we multiply the APEC 
model by the PLABS model, E~ a , with a = -0.03 and fit it 
to the XIS+PIN spectra. The best-fit parameters are obtained 
to be kT = 11.85 keV, Z = 0.30 solar, and K = 1.35 X 10 -2 
(X 2 /dof - 1343/1219). This indicates that the Suzaku tem- 
perature measurement has a systematic error at the level 
of ~ 1 keV, however, the discrepancy between Suzaku and 
Chandra becomes even larger in this case. Thus the uncer- 
tainty of the Suzaku effective area is unlikely to be the major 
source of the discrepancy. 
2. Chandra ACIS-I 

Based on a comparison of cl uster temperatures de rived from 
Chandra and XMM- jVewfow. ISnowden et all (120081) reported 
that the Chandra temperature tends to be by 1-2 keV higher 
than that of XMM-Newton particularly for the clusters with 
kT > 5 keV. Furthermore, the proceedings of the Chandra 
Calibration Workshop by David et al. Q indicated the same 
tendency as well as that the temperature obtained from the 
fit to the 2-6 keV continuum spectra is systematically higher 
by about 2 ke V than that inferred from the iron line ratios for 
hot clusters. Thus we consider that the above results suggest 
that the Chandra effective area has a systematic error of the 
level of AkT = +1 ~ +2 keV for hot (kT > 5 keV) clusters. 
For example, fitting the Chandra data with the APEC model 
multiplied by the a = -0.03 PLABS model leads to 
kT = 12.89 keV, Z = 0.40 solar, and K = 1.44 x 10~ 2 
(x 2 /dof = 569/426). Although the metal abundance is 
marginally higher, the temperature agrees with that of XIS. 
Then by simultaneously fitting the Chandra and XIS spectra 
under the APECxPLABS model and determined a and the 

relative normalization factor to be -0.017(-0.031 0.003) 

and 1.07(1.06 - 1.08), respectively. The best-fit APEC pa- 
rameters are kT = 12.92 keV, Z = 0.35 solar, and K 
1 .37 x 10~ 2 (y 2 /dof = 1890/1636). The result is plotted in 
Fig.lATTl 



Based on the above discussion, in order to take into account 
the difference of the calibration between the two satellites, we 
assume in §0 the Chandra normalization factor relative to XIS 
of 1 / 1 .07 = 0.93 and multiply the spectral model derived with 
Chandra by the PLABS model with a = 0.02(0 - 0.03). 



7 http://cxc.harvard.edu/ccw/proceedings/07_proc/presentations/david/ 



